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Summary. — The recent analysis of the cosmic microwave background data carried 
out by the WMAP team seems to show that the sum of the neutrino masses is 
< 0.7 eV. However, this result is not model-independent, depending on precise 
assumptions on the cosmological model. We study how this result is modified when 
the assumption of perfect lepton symmetry is dropped out. 

PACS 95.35. +d - Dark Matter. 
PACS 98.80.-k - Cosmology. 
PACS 14.60 -z - Leptons. 



1. Introduction 

The anisotropy spectrum of the cosmic microwave background is without doubt one of 
the most important observables quantity in cosmology. Recently, its precise measurement 
by the WMAP experiment, has lead to a precise determination of the values of the 
parameters describing our Universe Q . In particular, the WMAP team finds a very tight 
bound on the sum of neutrino mass, representing an improvement of roughly one order 
of magnitude with respect to the laboratory limit. It seems that finally the exciting 
possibility of performing precision tests in particle physics models using cosmological 
data, is within our reach. 

Unfortunately, parameter estimation is plagued with the problem of parameter de- 
generation. With this it is meant that different parameters have similar effect on the 
microwave background spectrum; the latter cannot be used, then, to distinguish between 
them. A well known example is the degeneracy between the spectral index n and the 
reionization redshift z nr . Put in another way, we could say that there is not a one-to-one 
relationship between the parameters and the spectrum. The simplest way to break such a 
degeneracy is to measure some other quantity; for example, the matter power spectrum. 

Another difficulty in dealing with the degeneration of parameters is due to the fact 
that we do not have an explicit relationship between the parameters and the spectrum. 
This means that probably there are some degeneracies we are not aware of, and that, even 
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for kwown degeneracies, we do not know exactly how much they influence our results. 
Thus we are forced to use numerical methods, that are in some cases very time expensive. 

Our work deals with the possibility to extract from the microvave background infor- 
mations on the degree of lepton asymmetry existing in our Universe. In particular, we 
study how the estimation of cosmological parameters is affected by the existence of such 
an asymmetry. The paper is structured as follows: in section 2, we review the basis of 
neutrino physics; in section 3, we briefly recall the basis of the microwave background 
physics ; in section 4, we describe the techniques for parameter estimation; in section 5 
we speak about the role of neutrinos in cosmology; in section 6 we present and discuss 
our results; finally in section 7 we draw our conclusions. 

2. Neutrino Physics 

In this section we review some basic facts of neutrino phyics, with particular regard 
to the formalism and to the phenomenology of neutrino mixing. 

Neutrino oscillations are due to the fact that the flavour eigenstates v 1 = {v ei i>^, v T ), 
defined on the basis of the interaction with the charged leptons (so that v e is the partner 
of the electron, and so on) do not coincide with the mass eigenstates (y = u\, V2, ^3), 
but are instead related by a unitarity transformation: 

(1) v' = Uv 

where U is the unitary 3x3 mixing matrix, parametrized in terms of three mixing angles 
(#12, #23, #13) and one phase (f>. Denoting mdi ag = Diag (mi, TO2, ^3), the effective mass 
matrix m v is given by: 

(2) m v = Um d - mg U T 
It is customary to label the mass eigenstates so that: 

(3) 0<Am| 1 <|Am^ 2 | 

where Am| = \m\f — \m\j. The frequency Wjj of neutrino oscillations is determined by 
their mass difference: 

I A-m 2 . I 

(4) Uij 



AE 



The two distinct frequencies measured in the experiments involving atmospheric and 
solar neutrinos correspond, in terms of the mass eigenstates, to: 

(5) Am 2 sun = Am 2 21 , Am^ tm = |Am§ 2 | 

At the present time, we have the following information about the mass square differ- 
ences 0: 

(6) Am 2 un - 0.7- 10" 4 eV 

Am 2 - 0.2 • 1CT 2 eV 



EFFECT OF NEUTRINO ASYMMETRY ON THE ESTIMATION OF COSMOLOGICAL PARAMETERS3 



and the following about the mixing angles 0: 

(7) 6»i2 ~ 30° - 35° 

023 ~ 45° 

013 < io° 

However, we are still missing some important informations: first of all, the absolute 
values of neutrino masses, and consequently the mass hierarchy. Given the experimental 
data, there are three possible scenarios: 

\mi \ ~ \m,2\ ~ |m.3| 3> \rrii — rrij\ Degenerate 

(8) \ m i \ ~ l TO 2| | ttt-3 | Inverted hierarchy 
|m3 1 3> |mi| ~ |ma| Normal hierarchy 

Secondly, we need a more precise determination of the three mixing angles, in partic- 
ular #i3 for which only an upper bound is known. Finally, it is still unknown if neutrinos 
are Dirac or Majorana particles. If the latter one is the case, two additional parameters, 
namely the relative phases <f>\, <f>2 between the Majorana mass terms mi, to.2 and 1713, 
are needed. 

The tritium (3 decay experiments give the following bound for m e Q : 

(9) m e <2.2eV 

It should be noted that, since u e is not a mass eigenstate, its mass can only be defined 
as an expectation value: 



(10) m 2 e =J2\Ue 



|2 2 
TO,- 



Another important, complementary experimental approach is the search for the neu- 
trinoless double (3 decay (in brief, 0v(3f3 decay) . This process is sensitive to the "effective" 
neutrino mass TO ee , namely the 11 entry of the mass matrix to„: 



(11) 



3 
i=l 



E^e 2 



The evidence for §v(3(3 decay, if found, would imply among others that lepton number is 
violated and that neutrinos are Majorana particles. At the present only an upper bound 
on I TO ee I is known: 

(12) |TO ee |<0.5eV 

where we have quoted the more conservative limit. This bound cannot be directly trans- 
lated into a bound on m e , since the relation between this and TO ee is not unique, as it 
can be seen by comparing eqns. IjlOJl and |H) , but instead depends on the values of the 
mixing matrix elements. This fact can be used to distinguish between different models. 
In particular, if TO ee is found to be close to the upper bound (|12J) . as it has been recently 
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claimed, the inverse and normal hierarchy scenarios would be strongly disfavored with 
respect to the degenerate one. 

Finally, very recently the analysis of the cosmic microwave background anisotropy 
spectrum as observed by the Wilkinson Microwave Anisotropy Probe (WMAP), combined 
with the data on the galaxy power spectrum, has lead to the following bound on the sum 
of the three mass eigenvalues, roughly an order of magnitude more restrictive than the 
laboratory limit Q: 



It should be noted however that, since many parameters contribute in giving to the cos- 
mic background radiation and galaxy spectra their actual shape, this bound is model 
dependent, relying on several assumptions on the values of the other cosmological pa- 
rameters, such as the total density f2. In effect, more conservative analyses lead to a 
somewhat higher limit, of the order of 1 eV. In the following we shall discuss in more 
detail the issue of the stability of the cosmological bound on neutrino masses when some 
assumptions, mainly the neutrino-antineutrino symmetry, are relaxed. 

3. Cosmic Microwave Background Radiation 

The cosmic microwave background (hence thereafter simply CMB), namely the 2.7 K 
black body radiation that fills our Universe, was originated when baryons and photons 
decoupled and the latter began to propagate freely toward us. For this reason the CMB 
encodes a wealth of information about the Universe as it was at the epoch of decoupling, 
roughly 300.000 years after the Big Bang, corresponding to z^ cc ~ 1100. 

The first striking feature of the CMB radiation (apart from it being a nearly perfect ( 1 ) 
black body) is its extremely high degree of isotropy: the temperature varies by only a 
few parts in 10 5 all over the sky. This means that the early Universe was indeed a very 
homogeneous and isotropic place, being well approximated by a Friedmann-Robertson- 
Walker (FRW) model. 

However, much attention has to be paid to the small anisotropies in the CMB spec- 
trum, since they are strictly related to the density perturbations at the epoch of de- 
coupling that, growing, have originated the structures we observe today. The anisotropy 
spectrum can then be used to put very precise constraints on the overall properties of the 
Universe, i.e., to the cosmological parameters, since, as already noticed, the Universe at 
Zdoc is a conceptually simple place, well described using the linear theory of perturbations 
on an homogeneous and isotropic background space-time, and the physical processes rel- 
evant at that time (mainly the Compton scattering that keeps matter and radiation at 
equilibrium) are very well understood. In this section we describe how the anisotropies 
are generated; in the following we will explain how they can be used to estimate the 
cosmological parameters. 

3T. Causes of temperature fluctuations. - The temperature fluctuations that we ob- 
serve when looking at the CMB in different points in the sky, are usually classified 



( 1 ) It is worth stressing the main limitation in detecting the deviation of the CMB spectrum 
from a black body one is our capacity of building a reference black body. 
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according to their origin, in the following way. The primary anisotropics arc produced 
at the time of the decoupling; we will be mainly concerned with them. The secondary 
anisotropics are produced by processes occurring in the way of the photons from the 
last scattering to us. One example is the integrated Sachs- Wolfe effect (ISW), occurring 
when the gravitational potential in which the photon moves is varying with time; this is 
the case in a A dominated Universe. Finally, the tertiary anisotropics are due to the fact 
that the CMB is not the only source of microwave radiation in the sky, so that its signal 
is contaminated by several foreground emissions. The most important are the emission 
from the galactic dust, the synchrotron emission from the electrons accelerated in mag- 
netic fields, the thermal bremsstrahlung emission from hot electrons in the interstellar 
gas, and finally the emission from extragalactic point sources (AGN, blazars, quasars, 
...). Although the tertiary anisotropies are not anisotropics intrinsic to the CMB, so 
that their contribution has to be subtracted when looking for the cosmological signal, 
nevertheless they contain valuable information on several astrophysical phenomena. 

3 '2. Generation of primary anisotropies . - The primary anisotropies are produced by 
three basic mechanisms. The first one is directly related to the fluctuations in the matter 
density. Since before decoupling photons and baryons are tightly coupled, a fractional 
change in the density of the latter produces a change in the density of the former, and 
this in turn produces a change in the temperature of the radiation: a denser region will 
look hotter, while a less dense one will look colder. 

The second mechanism is the Doppler shift of the radiation. The peculiar velocity 
field of the baryons in the plasma is perturbed in a similar way to the density field; 
photons scattered by baryons that are moving toward us will look hotter, and vice versa. 

The third mechanism is the Sachs- Wolfe (SW) effect. Photons coming from overdense 
region will be redshifted due to the gravitational rcdshift, and then look colder: those 
coming from undcrdense region will be blueshifted and look hotter. 

Taking into account all these contributions, the equation for the fractional tempera- 
ture fluctuation AT/Tq at a direction n is: 

AT 1 1 

(14) — = -S 1 -h-v+- ( f> 

where To is the mean temperature, 5 1 is the fractional fluctuation in the photon energy 
density (the factor 1/4 coming from the fact that for a black body p 1 oc T 4 ), v is the 
velocity of baryons, <j> is the gravitational potential, and all quantities are evaluated at 
the time of decoupling. It is the combination of these three terms that determines the 
main features of the field of temperature fluctuations; since they are model dependent, 
the pattern of anisotropies can be used to distinguish between different models. 

3'3. The Power Spectrum. - The quantity used to describe the statistics of tempera- 
ture fluctuations is the power spectrum. The temperature field in the sky is expanded in 
spherical harmonics: 

(15) ^P=£ E a lm Y lm (n) 

1=1 m— — l 
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where 

(16) aim = y J dfiT (^ Y i™(™) 

are the multipole coefficients. The temperature angular power spectrum Ci is then given 
by: 

(17) («L«rm') = Ci5w5 mm i 

Where the angle brackets (...) denote an average over a statistical ensemble. Unfortu- 
nately, our sky is just one realization of such ensemble; however, owing to the fact that 
the spectrum does not depend on m, we can construct an unbiased estimator of C; in 
this way: 



(18) Q = ^— ]T at m a h . 



A given cosmological model will predict the values of the C; . This can then be used as 
a test for cosmological models. 

4. Parameter estimation techniques 

The problem of parameter estimation is often stated stated as follows: given a set 
of experimental data, we want to find the set of parameters that provides the best fit 
to the data themselves. In other words, we want to find the set of parameters that 
has the maximum probability to produce the observed data. However this is just part 
of the story, since the probability that a given model produces the data can be put in 
relation with the inverse probability, namely the probability that the data are produced 
by a given model, only once we assign some a priori probability do the different models 
(Bayes' Theorem). 

Thus, defining a maximum likelihood function C(d\ 9), being equal, apart from a nor- 
malization constant, to the probability that the model defined by the vector of parame- 
ters ( 2 ) 9 produces the data d, Bayes' theorem states that: 

(19) V{9\ d) oc £(d| 9) ■ V{9) 

where V(9\ d) is the probability that the data d come from the model 9, while V(9) is the 
a priori probability, usually called prior, assigned to the model 9. If we are interested only 
to one parameter di and not to the whole vector 9, we have to perform a marginalization 
of the distribution function over the remaining parameters: 



(20) P(6» l |d)oc J P{9\d)d6 1 ...de i - 1 de, 



The maximum likelihood principle states that the best estimator 9 of the set of pa- 
rameters 9 is the value for which V{9\&) assumes its maximum value. The problem 



( 2 ) In the following we shall simply say that 6 is "the model". 
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of parameter estimation is then reduced, at least in principle, to the knowledge of the 
function C and to the search of the maximum of the function defined in (|19|l with respect 
to the variable 6, being d known and fixed. In practice, however, such a search is not 
a simple task, in particular when one has to deal with a high dimensional parameter 
space. In the following we shall describe two techniques for the search of the maximum 
of the maximum likelihood function, supposing to be always able to compute £(d\ 0) for 
a given value of 0. 

4'1. 'Brute force' approach. - In this approach, the parameter space is sampled using 
a grid (see @ for an example of an application to the CMB data) . Let D be the number 
of dimensions in the parameter space (equal to the number of parameters) and let us 
consider, for the sake of simplicity, the same number n of intervals along each dimension: 
the total number of points in the grid is then n D . The function £ is then computed at 
every point in the grid and is multiplied by the priors. Let us assume that the priors are 
uniform, so that we can 'confuse' £(d| 0) e V(0\ d); the following is easy generalized to 
the case of non-uniform priors. Now, to obtain an estimation of the value of one of the 
parameters, we have to marginalize over the others, as said before. A direct integration 
is, however, computationally very expensive, and in many analyses it is substituted by a 
maximization over the remaining parameters. Let us consider for example a parameter 
to = Oi (i = 1, D); let ojj be the values of ui sampled over the grid (j — 1, ...,n). The 
procedure is then the following (for simplicity purposes we omit the dependence from 
the data d in the argument of the distribution functions): 

1. We fix the value of lo to a given u; 

2. We search for the maximum of C(0) for to = lu and assume that: 

V(uj = ui) oc max {£(#)} 

UJ— w 

3. We repeat the procedure for all the values ojj {j — 1, ...,n) and thus obtain a n- 
point sampling of ViuS)\ 

4. We perform an interpolation to obtain a likelihood curve for lu; 

5. We repeat the procedure for the others parameters. 

Now we have a probability distribution function for every parameter, and the problem 
is reduced to the one of finding the maxima of n functions of one variable, surely easier 
than its multi-dimensional analogous. 

The theoretical reason for the approximation described at the point 2) lies in the fact 
that the result of the marginalization (|20|) and that of the maximization are exactly the 
same if the maximum likelihood function is a multivariate gaussian. 
The main shortcoming of this approach is the fact that it requires large computational 
times. In the case of the CMB spectrum, the computation of the maximum likelihood 
function requires from 10 seconds to slightly more than a minute, on a typical workstation. 
If we consider a typical parameter space, with 7 dimensions, sampled in 10 point in 
each dimension, we obtain a grid made of 10 points. Assuming we need 30 seconds 
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to compute every model, the requested time is nearly 10 years. Using more powerful 
machines can reduce the computational times by some orders of magnitude, but the 
effort stays prohibitive for 10, or more, dimensional parameter spaces. 

4'2. Markov Chain Monte Carlo method. - The Markov Chain Monte Carlo method 
(simply MCMC) makes easier the search for the maximum of the likelihood function, 
using a 'clever' algorithm to sample the parameter space. It allows even, through the 
same algorithm, to compute with a lesser effort the marginalization integrals, avoiding 
in this way to approximate them with a maximization. Finally, it makes possible to 
explicitly compute the values of the normalization constants that we have overlooked in 
the preceding section. We shall briefly illustrate how this results are achieved; for further 
explanation see [y] . 

The idea, upon which the MCMC method is based, is that of building up a Markov 
chain whose target distribution is the one we want to sample (in our case, the V(9)). 
This means that, once the Markov chain has converged, it becomes a random number 
generator distributed according to the V{0). In this way it is possible to effectively sample 
the V(0); the Markov chain will 'visit' more often the points with higher probability, in a 
way proportional to the value of V in the point. The building of a Markov chain having 
as its target distribution the one we are interested to is made using a Metropolis-Hastings 
algorithm. 

It is the clear that the main advantage of the MCMC method is that the research 
algorithm does not waste time in zones of low probability. This implies that the compu- 
tational times scale with the number of parameters in a slower than exponential way. 

5. Neutrinos and Cosmology 

Our Universe is filled with a background of relic, non-interacting neutrinos, whose 
number density is of the order of 100 particles/cm 3 per species. Since, as pointed out 
by the recent experimental evidences, neutrinos are massive, they can give a non negli- 
gible contribution to the total energy density of the Universe. This contribution can be 
expressed, for a single species of mass m, as follows 0|: 



where h is the Hubble constant in units of 100 km sec 1 Mpc 1 . The function A(£) is 
defined as: 



where r/(n) is the Riemann 77 function of index n, and £ is a dimensionless parameter 
representing the chemical potential of the neutrino gas, normalized at the decoupling 
temperature. A non null chemical potential implies a lepton asymmetry, so we can think 
of it as quantity parameterizing such an asymmetry. In the case of no lepton asymmetry, 
equation 1)21(1 reduces to the well known Gerstein-Zel'dovich formula. 

The existence of a background of light (m<l Mev) neutrinos acts in a different way 
on the observable cosmological quantities, that can then be used to extract informations 
on the properties of neutrinos Between others, one can get informations on the mass, 



(21) 



n u h 2 = 



rn 



A(0 



93 cV 



(22) 
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on the lepton asymmetry, on the number of families, other than on the existence of sterile 
neutrinos, of radiative decays, and on the neutrino mean life. 

In particular, the CMB spectrum is affected by the quantity /„ = Sl v /Q m , namely 
the ratio between the neutrino density and the total matter density. In effect, since 
neutrinos are ultra-relativistic at recombination, their presence shifts the time at which 
matter-radiation equality occurs, a parameter whose value affects the CMB spectrum. 
As we can see from eqn (|21|) . f v depends on the mass, the chemical potential and the 
number of families. 

The matter power spectrum is itself affected by j Vl since the small-scale fluctuations 
are suppressed in a degree proportional to /„, due to the high velocity dispersion (free 
streaming) of the neutrino gas. 

Finally, the abundance of 4 He produced during the primordial nucleosyntheys is af- 
fected both by the the total density of neutrinos, since this influences the expansion ve- 
locity and consequently the reactions freeze-out, and by the neutrino-antineutrino asym- 
metry in the electronic sector. Combining the nucleosynthesis data together with the 
CMB data, the following bounds are obtained [|J: 

(23) l^rl < 2.6 

(24) 1^1 < 0.3 

On the other hand, if the Large Mixing Angle solution for the mixing matrix is true, the 
oscillations equalize the chemical potentials and then 

(25) |^,-r,e| < 0.07 

In this way, by obtaining and independent estimation of £, it would be in principle 
possible to falsify the LMA solution. 

About the mass, the best estimation is the obtained using together the CMB data of 
the satellite Wilkinson Microwave Anisotropy Probe (WMAP), of the Cosmic Background 
Imager (CBI) e of the Arcminute Cosmology Bolometer Array Receiver (ACBAR), the 
data on the large scale structure of the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS) and the data on the matter power spectrum at z ~ 3 obtained by the observations 
of Lya forest. The resulting upper bound, obtained in the hypothesis of three species 
with the same mass and null chemical potential, is (as already quoted in section [2|l }]'■ 

(26) m v < 0.23eV 

As we have already noticed, this represents a improvement of one order of magnitude with 
respect to the laboratory limits. However, some authors have criticized this result [TlHl2] |. 
on the basis of the fact that the interpretation of the Lya data is still controversial, and 
of the fact that there are some degeneracy effects in the paremeter space. In particular, 
in 0| the influence of the priors on fi m and on the bias parameter b has been discussed 
, and it has been shown that the WMAP and 2dFGRS data are well fitted by a model 
having Q u = 0.2, Q m = 1 e fl\ = 0, altough this would require a very low value for the 
Hubble constant, h < 0.5. In a similar way, in the authors stress the degeneracy 
existing between the neutrino mass, the Hubble constant and the bias parameter. The 
latter is a phcnomcnological parameter that is introduced to take into account the fact 
that the galaxy distribution could not trace exactly the actual matter distribution, and 
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then be biased. The main advantage in combining the CMB data with the data of the 
large scale structures, is that the former can be used to fix the amplitude of fluctuation, 
and then to put some bound on the bias. However, since as said before the effect of 
neutrinos is the suppression of the power on small-scales, the two effects can be in some 
limit confused. 

Finally, the WMAP analysis assumes that the chemical potential of all species is equal 
to zero, and then that a perfect symmetry exists in all leptonic sectors. From the the 
theoretical point of view, there is no reason to suppose this. 

6. Effect of a lepton asymmetry on the parameter estimation 

Nearly all the analyses of the cosmic background spectrum assume, often implicitly, 
a perfect neutrino-antineutrino symmetry, and the a null chemical potential. Some ex- 
ceptions exist, regarding the use of the CMB together with the abundance of 4 He to put 
a bound on the chemical potential; however it is assumed that the neutrino has a null 
mass [jj ^3 . In this way the effect of the chemical potential can be paremeterized in 
a simple way, introducing an effective number of neutrino families N Vte g. Furthermore, 
just two parameters are used: Qf,h 2 e N Vje g. In 0], the effect on the power spectrum 
of a non null neutrino mass and chemical potential has been studied, but the authors do 
not make a statistical analysis over the CMB data. 

It is then worth studying the effect that a prior on £ can have on the determination 
of the other parameters, and trying to understand if there exists a degeneration between 
it and the other parameters. For this, we have performed a preliminar analysis of the 
BOOMERanG 98 (B98) data, using the techique described in the section ITT1 In the 
following we shall describe it in more detail. 

The first step is the choice of the parameter space. We have restricted our attention 
to flat models, with baryons, cold dark matter, cosmological constant and three neutrino 
species with the same mass and chemical potential. We have 6 parameters: the physical 
baryon density f2&/i 2 , the density parameter associated to the cosmological constant JIa, 
the scalar spectral index n, the Hubble constant h, the neutrino mass m„, the neutrino 
chemical potential £„. The values are listed in table |U 



Parameter Values 



n b h 2 

h 



{0.004, 0.009, 0.02, 0.045, 0.1} 
{0, 0.25, 0.5, 0.75, 1} 
{0.6, 0.7, 0.8, 0.9, 1.0} 
{0.90, 0.95, 1.0, 1.05, 1.10} 
{0.02, 0.033, 0.056, 0.093, 0.15, 0.26, 0.43, 0.72, 1.20, 2} eV 
{0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0} 



n 



Table I. - Parameter values 



The second step requires to compute the power spectrum for every combination of 
the 6 parameters. The total number of possible combination is about 70,000, but some 
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of them correspond to non-flat models. Eliminating these, we are left with about 50,000 
models. The spectra are computed using the CMBFast code, modified according to 01 
to take the chemical potential into account. 

The third step requires to process the theoretical spectra to obtain quantities compa- 
rable to the B98 data. For this we have to reduce the Ci to band power, averaging over 
near, correlated multipoles through a window function. In this way we obtain uncorre- 
cted points. The quantities obtained are called 

The fourth step is the computation of the likelihood of each model. We use the offset 
lognormal approximation introduced in |l5| : 



and Xb, o~b, Gbb' are known quantities. Finally, the last step is the marginalization, using 
the maximization method described in section im First we have considered only models 
having £ = 0, imposing then a prior on the chemical potentail; then we have considered 
the whole space. In this manner we obtain two likelihood curves for each parameter; we 
report them in in figs. 1 and 2. In fig. 1, we have just plotted the sampled points, and 
joined them by lines. In fig. 2 we show a tentative interpolation. 

The comparison between the two curves provides useful information on the effect of 
lepton asymmetry on the cosmological parameter estimation. We expect to have, in the 
case £ = 0, an agreement with WMAP's best fit values (but please note that we did not 
include the optical depth r in our parameters space) . Now let us discuss separately each 
parameter. 

The physical baryon density u>b = f^/i 2 shows a maximum at about 0.03-0.04, both 
for £ = and for freely vaying £. In this last case, however, the low baryon density 
models are less disfavoured than in the former case. WMAP's best fit value is 0.024. 

The vacuum energy density parameter S1a shows, in the case £ = 0, a clear peak 
at about 0.7. When we allow £ to vary, the curve get flat between and 0.7. Such a 
behaviour seems to point to a degeneracy between £ e In 16], a similar effect has 
been observed, with respect to the data on the large scale structures. It seems then worth 
looking deeper into this problem. WMAP's best-fit value is 0.7. 

The dimensionless Hubble constant h shows, when £ = 0, a maximum between 0.6 
and 0.8. When £ can vary, the maximum is shifted toward the left. WMAP's best-fit 
value is 0.72. 

For the scalar spectral index n, low [n < 0.95) values are favoured when £ = 0; when 
£ is freely varying, higher (n > 1.05) values seem to be preferred. WMAP's best fit value 



For the neutrino mass m„, low (m v < 0.5 eV) values are favoured when £ = 0; when 
£ is allowed to vary, a peak appears at 1 eV. The WMAP team put an upper bound to 
the neutrino mass of m„ < 0.23 eV. 

The dimensionless chemical potential £ shows two maxima: the first at 1.5, the second 
at 4. We note that the latter is not in agreement with the upper limit of 2.6 imposed by 
the primordial nucleosynthesis. 



(27) 




where 



(28) 



(b = ln(C& + Xb) 



is 0.99. 
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£2 a 



0.01 0.1 0.2 0.4 0.6 0.8 1 0.6 0.7 0.8 0.9 1 



n m v § 




0.9 0.95 1 1.05 1.1 0.1 1 1 2 3 4 5 



Figure 1. - Sampling of the likelihood curves for the six parameters we have considered. Upper 
row, from left to right: the physical baryon density u>b, the vacuum energy density Qa, the 
Hubble constant h. Lower row, from left to right: the scalar spectral index n, the neutrino mass 
m v , the chemical potential £. The thick lines correspond to the case with £ = 0, the thin lines 
correspond to the absence of assumptions on f. 



Finally, in fig [3] we show the maximum likelihood surface for the pair {m„, £„}. We 
see three maxima, one at {0.5 eV, 1.5}, one at {1.5 eV, 1.5} and one at {1 ev, 5}. This 
last one would be excluded on the basis of the upper bound cited above. 

7. Conclusions 

The analysis described in the preceding section points to the fact that the assumptions 
usually made on the lepton asymmetry do really influence the estimation of cosmological 
parameters. In particular, the parameters that are more affected are the cosmological 
constant, the neutrino mass and the spectral index. Values of the cosmological constant 
between and 0.7 cannot be distinguished, having nearly the same probability. The 
neutrino mass is peaked at 1 eV, this being not in agreement with the WMAP results. 
Finally, higher values of the spectral index are favoured. We think then that is worth 
studying in more depth this topic. It is without doubt necessary to apply our method to 
the WMAP data, using a larger parameter space, to be explored using MCMC methods. 
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